Radar reflection off extensive air showers 
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Abstract. We investigate the possibility of detecting extensive air showers by the radar 
technique. Considering a bistatic radar system and different shower geometries, we sim- 
ulate reflection of radio waves off the static plasma produced by the shower in the air. 
Using the Thomson cross-section for radio wave reflection, we obtain the time evolution 
of the signal received by the antennas. The frequency upshift of the radar echo and the 
r-j \ power received are studied to verify the feasibility of the radar detection technique. 

1 Introduction 

Traditional techniques of extensive air shower (EAS) detection include counting the shower particles 
at the ground level and measuring fluorescence light from the excited nitrogen in the atmosphere. 
Detecting radio emission from the shower particles is a promising new technique that is currently in 
development. An alternative method is the radar technique, where a ground-based radio transmitter 
illuminates the ionization trail left behind the shower front and another ground-based antenna receives 
the scattered radio signal with its upshifted frequency. 

The concept of implementing a radar for cosmic ray detection was first introduced about 70 years 
ago by Blackett and Lovell [1|. However due to the negative results of dedicated experiments, this 
concept was forgotten for several decades. In recent years renewed attention has been attracted to this 
topic, see e.g. B2I3I4I5I . Experimental efforts to detect EAS using the radar technique were made by 
several groups, e.g. @, the LAAS group Q, 0, the MARIACHI experiment [9 1, or the TARA project 
IflOll . however no detection was confirmed so far. 

In this paper we investigate the possibility of detecting EAS by the bistatic radar system. Es- 
pecially, we are interested in verifying the feasibility of supplementing the CROME (Cosmic-Ray 
Observation via Microwave Emission) detector 11 1 II 21 with a radio transmitter and using its existing 
antennas as receivers for the radar echo. The CROME antenna setup consists of several microwave re- 
ceivers for the frequency ranges 1 .2 - 1 .7 (L band), 3.4 - 4.2 (extended C band), and 10.7 - 1 1 .7 GHz 
(low Ku band). Since we expect large frequency upshifts of the radar echo and the CROME antennas 
are fine tuned to the GHz frequency band, we can narrow down our interest to the frequency range of 
1 - 100 MHz of the emitted radio waves. 

The concept of detecting EAS using radar technique is based on the principle of scattering radio 
waves off the static plasma produced in the atmosphere by the energetic particles of the shower. The 
locally produced plasma decays exponentially with the lifetime t. According to Vidmar iTPJl . for the 
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Fig. 1. A schematic diagram representing a bistatic radar system and reflection from the static plasma disk pro- 
duced by the EAS in the atmosphere. See the text for a detailed explanation. 



plasma densities relevant for EAS and at low altitudes, the three-body attachment to oxygen dominates 
the deionization process as it depends quadratically on the oxygen density. This leads to the plasma 
lifetime of 10 ns at sea level and about 100 ns at altitude of 10 km. Since the received power of the 
radar echo is strongly diminished by the geometrical factor, the strongest signal will be obtained from 
the altitudes close to the ground level. Thus we can safely assume t = 10 ns. 

The ionization trail that results from meteor or lightning is traditionally divided into the underdense 
and overdense regions, depending on the local plasma frequency v p . If the electron density is high 
enough that the plasma frequency exceeds the radar frequency, then the region is overdense and the 
radio wave is reflected from its interface. In contrast, if the local plasma frequency is lower than the 
frequency of the incoming radio wave, then the region is underdense and the radio wave can penetrate 
the ionized region. In such a case the reflections are caused by the Thomson scattering of the radio 
wave on individual free electrons. 

If we consider the radar frequency in the range 1-100 MHz and shower energies up to 10 20 eV, 
then the maximum size of the overdense plasma region will be several meters. Now, if we calculate 
the attenuation length of the radio wave in the plasma with frequency v p = 100 MHz and v p — 1 
MHz one obtains ~300 m and ~3000 km, respectively. Those lengths are much larger than the size of 
the overdense plasma, thus the radio wave can easily penetrate the whole volume of the plasma disk 
produced by the EAS. 

Using the Thomson cross-section for radio wave reflection seems to be justified from the above 
arguments. An exact calculation of the total power received requires the integration of the contributions 
from each individual electron in the shower. The final result will then depend on the individual phase 
factors of the scattering electrons. 



2 Radar reflection 

A schematic diagram representing a bistatic radar system is shown in Figure Q] A ground-based radio 
transmitter (T) irradiates a disk-like shower. The front edge of the disk is an ionization front moving 
towards the ground with the speed of light and leaving created static plasma behind it. The radio 
signal is scattered by the free electrons in the ionization trail and subsequently received by the ground- 
based antenna (R). The geometry of the bistatic radar system is set up by the polar coordinates of the 
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Fig. 2. The time evolution of the region from which the scattered radio waves arrive simultaneously to the detector. 
The arrival time t = corresponds to the moment at which the shower hits the ground. The considered shower 
has energy of 10 18 eV and is heading towards the transmitter (x=0 m, y=0 m). The receiver (x=200 m, y=0 m) is 
placed 200 m away from the transmitter. The colors scale with the altitude. 



transmitter and the receiver, i.e. by the distances from the shower core to the transmitter (dj) and to 
the receiver (<1r) together with the angles <p t and <p r . 

Let us consider the radar reflection from the disk element with coordinates (r^, <p), altitude h, and 
electron density n e . Its contribution to the radar echo at the receiver at time t can be expressed by 

U rcv (t, s, r Lt <p) = U T s/Gre 90 — — sin a J — \ n e , ( 1 ) 

|r||r sc | V dQ 

where Gj is the transmitter antenna gain, Uj [V] is the magnitude of the transmitted field, k and k sc 
are the wave vectors of incoming and scattered radio wave, <^o is the initial phase of the emitted signal, 
dcrjIdQ is the Thomson cross-section, s is the projection of the distance between the shower core and 
the considered disk element on the shower axis, a is the inclination angle of the reflected radio wave, 
and n is the refractive index of the air. By analogy to the setup of the CROME detector, we assume 
that the receiver is oriented vertically upwards. 

The signal received by the antenna at a given time t is sum of the signals scattered at different times 
and from different parts of the plasma disk. These individual contributions interfere with each other 
and only an integral over the whole volume uol(t) from which they arrive simultaneously gives us the 
correct value. The total electric field strength of the radio wave at the receiver is given by 

Urcv(t)= J U rcu (t,s,r L ,<p)r L dr L d<pds. (2) 

J [)£>/(/) 

Note that time t is defined in such a way that t — coincides with the moment at which the shower 
hits the ground. 
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Fig. 3. Upshifts f r of the radio wave scattered on different parts of the shower disk, which is heading vertically 
towards the transmitter. The altitude of the disk element is given by h, whereas its distance to the receiver in the 
horizontal plane is d. 



Figure[2]shows the time evolution of the region of integration uol(t) (see equation (0) for a shower 
with an energy of 10 18 eV heading towards transmitter (x=0 m, y=0 m). The receiver (x=200 m, y=0 
m) is placed 200 m away from the transmitter. The colors scale with the altitude. 

As we can see, the volume over which we integrate the signal spans usually over the wide range 
of altitudes. This implies a high variation of the air density. Therefore, the refractive index and the 
Cherenkov angle can vary significantly over the considered region. This makes the whole analysis 
very complex. 



3 Frequency upshift 

Despite the fact that the radio wave is scattered on the static plasma, the region in which the plasma 
is created moves with relativistic velocity. Therefore, we can observe the Doppler shift of the received 
radar echo. Figure [3] shows the upshifts f r of the radio wave scattered on different parts of the shower 
disk, which is heading vertically towards the transmitter. The altitude of the disk element is given by 
h, whereas its distance to the receiver in the horizontal plane is d. 

The frequency upshift depends on the wave direction and the refractive index of air. It has the 
highest value for the case in which the viewing angle coincides with the Cherenkov cone. The typical 
f r is high enough to upshift the MHz signal into the GHz range. Note that f r only weakly depends on 
the distance from the shower core to the transmitter dj. 



4 Details of the simulation 

The main assumptions and approximations are the following: 

- The shower longitudinal profile is taken to be a fit to the Gaisser-Hillas function for proton showers. 

- The lateral distribution of the shower is given by the Gora function lfT4l . 

- Static plasma produced by the shower decays exponentially with the characteristic time r = 10 ns. 

- We consider radar reflections only from the part of the ionization trail which is up to 40 ns behind 
the shower front. In other words, the length of the plasma disk produced by the extensive air shower 
is constrained to 4r. This leads to less than 2% of the total plasma electrons not being accounted 
for. 

- We use the US Standard Atmosphere model. 

- We assume that the initial phase of the emitted signal is tpo = 0. Choosing different values of ipo 
can lead to slightly different values of the received power. 
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Fig. 4. The waveforms R(t) of the radar echoes for different frequencies v of the emitted radio wave and differ- 
ent shower core-receiver distances d R . In all cases a vertical shower with energy 10 18 eV heading towards the 
transmitter is considered. 



- The effective area of the receiver antenna is Ar-1 m 2 and the transmitter emits signal over the 
whole upper hemisphere (Gt - 2). 

- The receiver points vertically and its sensitivity is independent of the direction and frequency of 
the radar echo. The detector efficiency is 100%. 



5 Power received by the detector antenna 

The ratio of the 'instantaneous' power P«(f) received by the detector antenna to the power emitted by 
the transmitter Pj is given by 

P R (t)/P T =R 2 (t) (3) 

where 

*4^A ,4, 

V47T Ut 

The term R{t) is proportional to the electric field strength detected by the receiver. It can be used in the 
Fourier analysis to obtain the power spectrum of the recorded signal. The 'real' power received by the 
detector Pr can be obtained by averaging R 2 (t) (according to the time resolution of the detector). The 
ratio Pr/Pt ~ GtA r . 



6 Results 

Figure|4]shows the waveforms R(t) of the radar echoes for different frequencies v of the emitted radio 
wave and different shower core-receiver distances cIr. In all cases the receiver is outside the Cherenkov 
cone (which is the most common case) and vertical showers heading towards the transmitter are con- 
sidered. The time t = coincides with the moment at which the shower hits the ground. 

As we can see, the frequency of the received signal is higher than that of the emitted one. Moreover, 
the amplitude grows with time, which is the geometrical effect of the reflections from the lower parts 
of the atmosphere. 

For showers, where the receiver is inside the Cherenkov cone, the time sequence is reversed and 
the lower part of the shower is seen first. Therefore, the amplitude decreases with time. However, this 
is the case only for very inclined showers with cores up to several meters away from the receiver. 

The ratios of the power received by the detector to the emitted one (Pr/Pt) for vertical showers are 
given in tableQ] It is evident that the strength of the signal decreases with the frequency v. It is simply 
because the size of the region from which we get the coherent signal decreases with the wavelength 
and the negative interference cancels out the signal. 
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Table 1. The ratios of the power received by the detector to the emitted one (Pr/Pt) in dB for different frequencies 
v of the emitted radio wave, different transmitter-receiver distances d R and different shower energy E. In all cases 
a vertical shower heading towards the transmitter is considered. 
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7 Conclusions 

The typical ratio of the received to the emitted power is in the range of 1(T 10 - 1(T 13 , which makes 
the detection of the radar reflection plausible. The C band (~3 GHz) is particularly well-suited for 
this purpose because of very low noise (< 10 K) in this frequency range. The CROME antennas 
supplemented with commercial high power MHz transmitter would create a radar system dedicated 
to EAS detection. A note should be added, that the shown time traces are for an infinite bandwidth 
detector - a realistic detector would only be able to detect the signal in narrow frequency range. 
Moreover, it will only see the shower for a certain fraction of its development. 

The performed analysis will help to choose the frequency of the emitted wave to optimize the 
detection of the radar echo. The work is in progress. 
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